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o : 

D , A kinetic one-dimensional Ising model on a ring evolves according to a generalization of Glauber 

■ rates, such that spins at even (odd) lattice sites experience a temperature Te {To). Detailed balance 
is violated so that the spin chain settles into a non-equilibrium stationary state, characterized by 
multiple interactions of increasing range and spin order. We derive the equations of motion for arbi- 
trary correlation functions and solve them to obtain an exact representation of the steady state. Two 
nontrivial amplitudes reflect the sublattice symmetries; otherwise, correlations decay exponentially, 

■ modulo the periodicity of the ring. In the long chain limit, they factorize into products of two- 
O ' point functions, in precise analogy to the equilibrium Ising chain. The exact solution confirms the 

expectation, based on simulations and renormalization group arguments, that the long-time, long- 
distance behavior of this two-temperature model is Ising-like, in spite of the apparent complexity of 
the stationary distribution. 
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^ ■ I. INTRODUCTION 

^ ' 

^ I Exact solutions of simple model systems play a key role in statistical mechanics. They provide reliable quantitative 
information about regular as well as singular behavior and serve as proving grounds for approximation schemes. 
Moreover, they can set important pointers in areas where a general theoretical framework is still lacking. A prime 
example is the study of non-equilibrium stationary states (NESS). Similar to systems in thermal equilibrium, NESS 
are characterized by time-independent macroscopic observables; however, at present these can only be computed with 
explicit reference to the imposed dynamics: there is, as yet, no equivalent of Gibbs ensemble theory for far-from- 
equilibrium steady states. As a consequence, most progress to date is made by studying specific models. 
\^ I Since the Ising model [1] is one of the most intimately known interacting many-particle systems, many non- 
^sj . equilibrium models depart from it by imposing an external force which drives the system out of equilibrium. Examples 
T-H ■ include couplings to multiple heat baths [2-4] , competing Glauber and Kawasaki dynamics [5] , or a current-inducing 
, global bias [6]. In all of these models, reviewed in [7], the non-equilibrium perturbation violates the detailed balance 
symmetry of the equilibrium dynamics. On large time and length scales, one can probe how this affects universal 
behavior, or, adopting a more microscopic but no less fundamental perspective, one can ask how non-universal prop- 
erties, such as the exact configurational probabilities, are modified. One finds, generically, that non-equilibrium forces 
C ' have especially profound effects if (i) they couple to the bulk rather than the boundaries [8,9] of the system [10], 
. and (ii) the dynamics satisfies a conservation law [11-13]. For example, Ising lattice gases with conserved particle 
(~| ' number can be driven into two distinct non-equilibrium universality classes, depending on the symmetries of the 
O ■ external forces [6,11]. In contrast, Ising-like systems with non-conserved dynamics remain in the universality class of 
^ ' the equilibrium Ising model, even if the usual Z2-symmetry of the Ising model is broken [14]. Even then, however, 
J> ' the violation of detailed balance leads to fundamental changes in the configurational probabilities which must be 
! ■ identified and interpreted before we may hope to formalize our understanding of NESS. It is here that exact solutions 
r> [ are expected to be most helpful. 

. Motivated by these considerations, we recently [15] investigated a very simple non-equilibrium Ising-like model, 
■ " " ' namely, a one-dimensional mteracimgr spin chain with spin- flip dynamics coupled to iwo temperature baths. The rates 
are a simple but nontrivial generalization of the familiar Glauber [16] rates: spins at odd (even) sites are coupled to a 
temperature To (Te). In two dimensions, this model [3,4] exhibits an order-disorder phase transition which belongs to 
the Ising class, according to renormalization group arguments [14] and Monte Carlo simulations [4]. In one dimension, 
the two-spin correlation function can be calculated exactly [17]. Its only singularity lies at To = T^ = 0, so that 
the lower critical dimension is d = 1. Seeking an expression for the steady state, we solved the master equation 
perturbatively, in an expansion in the temperature difference of the two heat baths, up to and including second 
order terms [15]. To our surprise, the full stationary distribution turned out to be rather complex: At each order. 
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additional spin operators appear, characterized by longer spatial range interactions and higher-order spin products, 
and lower-order coupling constants acquire corrections. So, at first order, one encounters next-nearest neighbor pair 
interactions, while at second order, two new terms appear: a next-ncxt-nearest neighbor pair interaction, and a four- 
spin interaction spanning four nearest-neighbor sites. All of these are allowed by symmetry, and none allowed by 
symmetry are absent. Given this structure, one can, at least in principle, extrapolate to higher orders in perturbation 
theory. 

Of course, these findings immediately raise an obvious question: how does this relatively complicated stationary 
distribution generate long-wavelength behavior in the Ising universality class? In this article, we pursue an alternate 
route towards the answer: Instead of aiming for the stationary state directly, we seek its representation in terms of 
correlation functions - since the knowledge of all stationary correlation functions is equivalent to knowing the steady 
state. Starting from the master equation, we derive a hierarchy of equations of motion for correlation functions of 
arbitrary numbers of spins. In d = 1, this hierarchy is closed and soluble. Remarkably, we find that - in contrast 
to the apparent complexity of the stationary distribution itself - the correlation functions are very simple. Apart 
from two nontrivial amplitudes which reflect the temperature difference between the two sublattices, the structure 
of correlations is completely analogous to the equilibrium Ising model: specifically, for A'^ — » oo, arbitrary (even) m- 
point correlation functions factorize into a product of m/2 two-point correlations. Of course, all correlation functions 
involving an odd number of spins vanish by symmetry. 

While some exact analytic results for steady-state distributions are available, they are confined to three classes of 
systems: first, one-dimensional lattice gas models, restricted to excluded volume interactions, such as the asymmetric 
exclusion process and its relatives [18 20], second, very special one-dimensional spin systems whose master equations 
are solved by the Ising Boltzmann factor [6,21,22], and third, interacting systems in one or two dimensions on very 
small lattices so that the number of degrees of freedom remains manageable [23]. To the best of our knowledge, 
the work presented here is amongst the first complete solutions for non-equilibrium stationary states with nontrivial 
(nearest-neighbor) interactions and arbitrary number of degrees of freedom. 

The article is organized as follows. We first introduce our model and its master equation. Next, we derive the 
equations of motion for arbitrary correlation functions. Following a brief review of the solution [17] for the two-point 
correlation function, we show how four-point correlations can be factorized into two-point correlations in the long 
chain ( A ^ oo ) limit. We then postulate that all correlation functions factorize in this manner, and show that 
this factorization solves the equations of motion. In an appendix, we generalize this solution to finite systems with 
periodic boundary conditions. We conclude with some comments and open questions. 



Our model is defined on a one-dimensional ring, with an even number A^ of sites, and periodic boundary conditions. 
A spin variable, ct, = ±1, denotes the value of the spin at site i, and nearest-neighbor spins interact according to the 
usual Ising Hamiltonian, 



with an exchange coupling J. The dynamics is a non-equilibrium generalization of the usual Glauber model [16]: 
spins on even and odd lattice sites experience different temperatures, Tg and Tq. Specifically, a configuration {a} = 
{ci, (72, ... , CTjv} evolves into a new configuration by flipping a randomly selected spin crj with a rate [17] 



II. THE MODEL 




(1) 



Wi{(Ti 




ai{ai-i + CTj+i) , 



(2) 



where 



J 7e = tanh(2J/fcBTe) , i even 
y 7o = tanh(2J/A;Bro) , i odd 



(3) 



Thus, the full time-dependent configurational probability p({cr}; t) evolves according to a master equation. 



N 



dtP{{cr};t) = ^ ^-Wi{, 



-(Ji)p ({ct}; t) -h Wi{-ai (7, 



(4) 



i=i 



where {a^^'} differs from {a} by a flip of the i-th spin. A trivial time scale has been set to unity, and we use 
dimensionless units for inverse temperature, i.e., /3e = J/{kBTe), etc. 
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Our goal in the following is to find a representation for the stationary solution ofEqn (4), q{{a}) = lim^^oc i)- 
This limit is unique, since Eqn (4) is ergodic: every configuration {a} can be reached in finite time from every other 
configuration {a'} (unless Te = Tq = 0). For equal temperatures T = Tg = To, the steady state is just the (canonical) 
distribution for the Ising chain, 

qo{W}) = ^exp{-H/kBT) . (5) 

It is of course straightforward to compute arbitrary correlation functions for the equilibrium Ising chain. Since 
the Ising model is invariant under a global spin flip (Z2-symmetry), only correlations of even numbers of spins 

are nonzero, and are easily expressed in terms of the parameter uj = tanh( J/ZcsT). Of course, we may - and 
always will - order the arguments of an m-point correlation function (m even) without loss of generality, such that 
1 < ki < k2 < ... < km < N. Then, in the N ^ oo limit one finds easily: 

{ak^...ak^y = a;(*^2-fei)+(fe4-fe3)+-+(fcm-fcm-i) 

= Wki(^k2y {(^ksf^k^y ■■■ {(^km-i^k^y^ for m even and TV ^ oo , (6) 

i.e., a general correlation function factorizes into a product of two-point functions. 

We now turn to the non-equilibrium model, characterized by different temperatures, T,, ^ Tg. The associated 
stationary state violates detailed balance [15] and differs from the Boltzmann distribution, Eqn (5). The degree to 
which detailed balance is violated can be measured by the parameter d = (70 — 7e)/2. Similar to the Ising model, 
the stationary state is invariant under a global spin flip (Z2). As a consequence, all stationary correlations of an odd 
number of spins vanish identically, and again, only even correlations need to be discussed. We also note the symmetry 
under translations by n lattice sites, T„, combined with d — > —d if n is odd: 

q{{a};d)=q{{-a};d) (7) 
g(M;d)=9(T„{(7};(-l)"d) 

While a direct exact solution of the master equation has proven elusive, a perturbative calculation, in powers of d, 
shows [15] that the stationary distribution for this non-equilibrium model is rather complicated, with longer-range and 
higher-order spin operators appearing. Specifically, writing q{{cr}) = exp [1/ ({cr})], we find that the potential 
function V ({tr}), to second order in d, has the form 

^ ({o"}) = X] + '^'^ <^»cri-l-2 + 2 {dXf ^ [fTjCri+2 + (JiCFi+i(Ji+2cri+i - coth(2^) {aiai+1 + CTjCrj+s)] , 

i i i 

. _ (8) 

where tanh(2/3) = (7e -I- 70) /2 and A = — i sinh(4/3). An analysis of the structure of the perturbation series indicates 
that, at each order, additional spin operators appear, consisting of spatially longer-range interactions and higher-order 
spin products. It is therefore quite remarkable that the correlation functions turn out to be very similar to those of 
the Ising model, as we now proceed to show. 



III. EQUATIONS OF MOTION FOR ARBITRARY CORRELATION FUNCTIONS 



We begin by deriving the equations of motion for arbitrary time-dependent m-point spin correlation functions, 
{'^ki'^k2---'^km)v starting from the master equation: 

N ( 
i=l [{a} 

Here, the subscript (o)^ will be used to distinguish time-dependent averages from their stationary limits, (o) = 
limt^oo (°)f Due to the sum over all configurations, the {} bracket on the right hand side obviously vanishes for all 
sites i which do not belong to the set {ki,k2, ...,km}, and one finds easily that 

jn 

dt (o'feiO-fe2—0'fe„)t = -2^ (a'feiCrfe2...crfe„Wfei(o-fei -<^ki))t ■ 
1=1 



-Wi{(7i -(Ti)p {{a}; t) + Wi{-ai 
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Inserting Eqn (2) for ?«fe;((Tfe. ^^-fei) and taking the infinite time limit, we obtain the equations satisfied by the 
stationary correlation functions, linit^^ {'^ki'^k2---'^k,„)t = {'^ki^'k^-'-^km)'- 

m 

= -2m(<TfciC7fe2...C7fe^) + ^7fe, (iTfei..crfc,_i(c7fe,+i +crfc,_i)crfc,+i..crfc^) (9) 

i=l 

As mentioned before, all stationary correlations of odd numbers of spins arc trivially zero, and we need to focus 
only on the case of even m. Moreover, all stationary correlations are translationally invariant, modulo the sublattice 
structure. As in the Glauber model [24], the hierarchy of correlation functions is closed: the equation for m-point 
correlations does not involve any higher correlations, and it is homogeneous provided there are no nearest-neighbor 
pairs among the arguments. Otherwise, the m-point correlations couple to lower correlation functions, which appear 
as inhomogcncitics. In the following section, we motivate an ansatz for {(Tki<^k2---'^km) i then show that it satisfies 
Eqn (9) for any choice of {fci, fe, fcm}- 



IV. EXACT SOLUTIONS FOR STATIONARY CORRELATIONS. 



To establish several key relations, we briefly review the exact solution for the two-point functions {cri<Jj) [17]. For 
this case, Eqn (9) reads 

= -4 {(JiUj) + 7i ((cTj+i + Ui--i)uj) + 7j {ui{aj+i + aj-i)) . (10) 
For nearest-neighbor sites, e.g., j = i + this equation becomes inhomogeneous: 

= -4 {uiGi+i) + 7i {ui-iUi+i) + 7i+i {(Tiai+2) 7i+i (11) 

For simplicity, we restrict the discussion in this section to correlation functions in the thermodynamic limit N ^ 00. 

Thus we may label the lattice sites by the integers, i ^Ij. Correlations on finite periodic chains and some details for 
the N ^ 00 limit will be addressed in Appendix A. Following [17], it is easy to show that the solution is unique and 
takes the form (for i < j, without loss of generality) 



= y^AUju;'-' (12) 
with the spatial decay length controlled by the parameter 

The key difference to the equilibrium Ising model is the emergence of two amplitudes, to match the sublattice sym- 
metries, namely, 

= (7e +7o)/(27o) , if i even , 
Ao = (7e+7o)/(27e), if i odd ^'^^ 

For later reference, we note that these relations imply a further identity for products of pair correlations, namely, for 
a.nyi<j<j+l<l, 

= (7j -I- 7j+i) (fJi ai) - -fj {ai aj+i) {aj+i ai) - 7^+1 {(Ji cjj) {aj ai) . (15) 

This relation will be needed in the discussion of four-point functions. 

It is quite remarkable that the three quantities w, Ag, and Ag also determine all higher correlation functions, 
through a structure which is almost perfectly analogous to the equilibrium case. Deferring the case of finite A'' to 
Appendix A, we consider only an infinite chain here. Focusing on even correlations, we assert that arbitrary m-point 
correlation functions, with fci < ^2 < ... < km and m even, are given by 

{<^k,...<rkj = ^Ak,Au,...Ak^ ^(fe,-fe,)+(/c4-/c3)+...+(/c„-fc„-i) = ^^^^^^j ... {okrn-.<^kJ) ■ (16) 

i.e., higher correlation fimctions factorizc into two-spin correlations as in the equilibrium case. Eqn (6). However, two 
features distinguish these correlations from their equilibrium counterparts. First, the spatial dependence is controlled 
by a different parameter, oj. In analogy to u) for the equilibrium system, w defines an effective temperature for the 
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non-equilibrium system, via lu = tanh( J/fesTcff )• ToH' diverges with To or Tg but vanishes only if both To and go 
to zero. As a result, the correlation length ^, defined via Incj = — diverges only if both temperatures vanish. 
Second, and more importantly, we note the appearance of the even/odd amplitudes A^, Ag, reflecting the sublattice 
identities of the two spins. These amplitudes carry the primary information about the non-equilibrium nature of our 
dynamics. 

Before turning to a general proof of Eqn (16), it is instructive to confirm it explicitly for the four-point functions. 
In this case, factorization implies 

{ci-,ajak(Ti) = (ffjO-j) {(TkO-i) , (17) 

for i < j < k < I. When inserting this ansatz into the right hand side of Eqn (9), we need to distinguish whether j 
and k are nearest-neighbor sites or not. If fc > j -|- 1, we obtain 

-8 {aiajakCTi) + 7j ((cj+i + (Ti-\)ajak(Ji) + 7j (a-j(crj+i + aj-\)ak(Ti) 

+7fc {(JiCrj{(Tk+i + (Jfe_i)cr/) + 7; {(Jiajak{(Ji+i + cr/_i)) 
= [-4 {aiaj) + 7j {{ui+i + (Tj_i)(Jj) + 7^- ((Ti(crj+i + o-j-i))] (cfeCT;) (17a) 

+ [-4 {(Tk<^l) + 7fe ((o-fe+l + <^k-l)<^l) + ll ((7'jO-fe((7(+i + Cr;_i))] (cTiCTj) 

= , 

since the expressions in the square brackets vanish for the stationary two-point functions, by virtue of Eqns (10) or 
(11), depending on whether (i, j) or {k,l) are nearest neighbors or not. If k = j + 1, we need to consider a slightly 
different equation, namely 

—8 {aiGjaj+iai) + 7^ (((Jj+i + ai-i)a jU j+iui) + 7^- {ai{aj+i + a j-i)a j+iai) 
+7j+i {(^i(^j{(^j+2 + o-j)(T;) +7/ (<7i(^jCTj+i(o-/+i + (Ti-i)) 
= [-4 (cTjCTj) + 7j (((Tj+i + (Ti-i)aj) + 7j {aiaj-i)] {uj+iui) + {aiui) 

+ i'^t'^j) [-4 {cTj+icn) + Ij+i {o-j+20-i) + 7/ ((Tj+i((j;+i + cr/_i))] + jj+i {(Jiai) (17b) 

= 7j (o-jfTi) + 7j + i {(JiCFl) - 7j (o-ifTj+i) (o-j + iCTi) - 7j+l (fTjCTj) (fTjCTi) 

= . 

The last equality follows from Eqn (15). 

The proof is completed by induction. Assuming that the factorization has been proven for (m — 2)-point functions, 
it is suSicient to show that the ansatz 

(o-fci...o-fc„) = (o-fci---o-fe„_2) (o-fc„,_iO-fc„) (18) 

solves the equations of motion for the m-point functions. Again, we need to distinguish whether km-2 and km-i are 
nearest neighbors or not. If they are separated by more than one lattice spacing, i.e., if km-i > km-2 + 1, we have 

m 

-2m {(JkiCrk2---'^km) {'^ki'-'^ki-A'^ki + l + ffe._i)fTfe._^i..fTfc^) 

m-2 

i=l 

+ {uki-(^km-2) [-2 (0-fc„_il7fe„) +7fcm-l ((o-fe^-l+l + I7fem-l-l)0'fem) + 7fcm (c^/Sm- 1 (^^m + l + f^fen. -1 ) )] (l^a) 




-2(m- 2) (a-fei...(7fe„_2) + ^ 7fe; (o-fei ..crft._i (crfe.+i + crft._i)crft._^i ..crft^.^ ) 



Here, we have used the fact that both square brackets vanish, since they enclose the equations of motion for the two- 
and (m — 2)-point functions, respectively. Next, we consider the case fcm-i = km-2 + 1- Again, we add and subtract 
the terms missing from the full equations of motion of the lower correlations functions: 

m 

-2m{aki(Tk2--0-k,n) + ^^^ki {o'ki-<^ki-i{o'ki + l + <^ki-l)o'ki+i--<^k^) 
i=l 

= ~7fem-2 (^fei ••0'fe„_3(Tfe„_2 + l) (crfc„_iCrfc„) + 7fc„_2 (tJfci ■•CTfc^.jCTfc^) 

-7fc™-i (o'fci--0'fc„_2) (crfc„_i-icrfc„) + 7fc™-i (c^fci--ffe„_3afc„) 
= (crfci..crfc„_J {(7fc„-2 +7fe™-i) (o-fe„-3C^fc„) (18b) 

-7fe„-2 (o'fc„-30'fcm-2 + l) (o'fc„_iCrfc„) - 7fc™-i (o-fe^.aCT/c^.j) (crfc„_i -iCTfc^ ) } 

= . 
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To obtain the last two identities, we have factored out an (m — 4)-point correlation, and used Eqn (15). This concludes 
the proof. We emphasize again that the condition ^ 1 has been imposed for reasons of simplicity alone: expressions 
for arbitrary correlation functions on finite rings, while somewhat more cumbersome, are easily derived and also follow 
the Ising pattern, as we will show in Appendix A. 

Finally, we discuss one special case of the two temperature Glauber dynamics where the two-point correlations do 
not have the form of Eqn (12). This is the case, if one of the two parameters 7e, jo is zero, the other finite, i. e. the 
temperature of one sublattice is infinite. Without loss of generality we assume in the following 70 = 0. One can then 
easily show directly from the equation of motion (10) that the two-point function reads 



Hence, spins are only correlated over a distance of at most two lattice sites. Turning to the higher correlations, we 
observe that Eqn (15) also holds for these parameter values. Therefore, m-point correlations (m even) again factorize 
into two-point correlations. 



To summarize, we have found an exact solution for all stationary correlation functions of a one-dimensional non- 
equilibrium Ising spin chain with an even number A'' of sites. The system is globally coupled to two temperature 
baths: spins on odd (even) lattice sites experience a temperature To (Tg) and flip according to a generalization of the 
familiar Glauber rates. The presence of two different temperatures violates detailed balance and maintains a nontrivial 
non-equilibrium stationary state. The complete set of correlation functions provides us with the full, exact solution 
for this stc^ady state. This allows us to reconcile two potentially contradictory earlier findings: while simulations 
and renormalization group arguments indicate that the long-wavelength behavior of non-conserved two-temperature 
models should be Ising-like, a perturbative analysis of the stationary distribution q{{a}) showed that a zoo of long- 
range and multi-spin interactions are present in \nq ({cr}). Here, we have shown that the exact correlation functions 
for this model are remarkably Ising-like: they decay exponentially with a characteristic correlation length which 
diverges only if both both To and Tg vanish. Also, in the N ^ 00 limit, arbitrary (even) m-point correlations factorize 
into products of two-point correlations, following exactly the same scheme in both the Ising and the two-temperature 
chain. The only key difference is the appearance of two nontrivial amplitudes, Ag and Ag, which reflect the sublattice 
symmetry: for each spin on an odd (even) site, the correlation function carries a factor of \/ Ao {\/A^). Of course, 
stationary correlations of odd numbers of spins vanish identically, due to symmetry. In short, the correlations of the 
model are entirely consistent with its Ising-like long-distance properties. 

It is remarkable - and not at all immediately obvious - that these two amplitudes should be the only remnants of 
the zoo of interactions in q ({c}). There is, however, a very simple and elegant representation [25] of the stationary 
state in terms of an extended Ising model, consisting of 2A^ spins on a comb-like d = 1 lattice with Hamiltonian 
Haux = -JJ2iSiSi+i- Jo J2odd i^i^i- Z^even , « i^-j • Here, the Spins { Sj } form an auxihary set which must be traced 
out in order to obtain stationary observables associated with the original variables {(Ji). If the interactions J, Jo and J^ 
arc tuned appropriately in the complex plane, the cr-correlations of TLaux arc identical to those of the two-temperature 
model. Further, the exact stationary state of the two-temperature model follows as q{{(y}) (x Trj^.j exp[— Houx]- 
Details and generalizations will be presented elsewhere [26] . 

It would of course be interesting to investigate other non-equilibrium versions of Glauber dynamics. Will the 
correlations still have a rather simple structure if a spin chain is coupled, in a translation-invariant manner, to n > 2 
different temperatures? Will they still factorize, provided a sufficiently large number of amplitudes is introduced? 
Clearly, we still need to cover lots of ground before even such simple non-equilibrium systems are fully understood. 
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In this Appendix, we address the correlation functions for our model on a finite periodic chain (a ring) of N sites, 
with i = 1,2, ...^N . To set the scene, we first review the equilibrium case, with uniform temperature T. The correlation 




(19) 



V. CONCLUSIONS 



APPENDIX A: STATIONARY CORRELATIONS IN FINITE CHAINS 
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functions of the Ising model on a ring are well known. Using an ordered set of arguments, 1 < fci < ^2 < ••• < km < N 
without loss of generality, one finds: 

(0-fei...(Tfe^)*^' = for odd m , 

{(Tkr-CThJ^" =_JL^{cj('=^-'=i)+('=*-'=3)+---+(fe™-fe™-i)+^^-(fe2-fei)-(fc4-fc3)-...-(^^ for even m . (Al) 

with u) = tanh( J/fcsT) < 1. The periodicity of the ring implies that {aki---crk^y is invariant under — > u)^~'^ for 
any integer < n < A^. For oo, Eqn (Al) obviously reduces to Eqn (6). 

Turning to the non-equilibrium case with two temperatures, we first assert that the stationary two-point correlation 
function, with I <i < j <N,]s given by 

= zR^^ H"' + ^'^"^'"^^] • 

where the parameter u) and the amplitudes Ai e {A^, Ao) are simply those of Eqns (13) and (14), and the normal- 
ization factor Z{ijj) is given by Z{uj) = 1-1- lo^ . Clearly, (ui aj) is periodic with period A''. For easy reference, we fill 
in some details here which are also relevant to the N oo case. 

Eqn (10) is essentially a second order difference equation in the variable j — i- To find a unique solution, we need two 
boundary conditions: one of these is Eqn (11) for j — i = 1 and the other is the requirement that (cTj aj) be periodic 
with period (or, for A'' oo, that {(Ji<Tj) vanish for large separations). We therefore expect a superposition of 
two linearly independent solutions whoso coefficients can then bo dotorniincd. Away from the boundaries, the ansatz 
(cTjCTj) = ^/AiAjU)^~^ reduces Eqn (10) to three relations, for the combinations i,j odd/odd, even/even, and odd/even 
respectively: 



= -4Ao + 2-fo\/Ae Ao (uj + uj-^) 
= -4Ae + 2je\/Ae Ao (uJ + UJ-^) 
= -4^/Ae Ao + (7oAe + 7eAo) (w + UJ-^) 

We note that the third equation is simply a linear combination of the first and second, expressing the fact that there 
are only two amplitudes, and not three, as one might have assumed (cf. [17]) based on the three types (odd/odd, 
even/even, odd/even) of pair correlations. Moreover, the symmetry under oj ^ uj~^ implies that, for each solution w, 
co~^ is also a solution. Proceeding to solve the system, we find 7eAo = joAg and 



W + W ^ = ^ LJ± = — ^= (l ± -7e7o) 



with a;_|_ = 1/ iO- . The two roots lo± of this quadratic equation provide us with the two anticipated linearly independent 
solutions. However, to satisfy the inhomogeneous Eqn (11) for both finite and infinite A'', we have to take a convex 
combination of the two solutions 

{uiOj) = Q!-y/A~A~a;i"* + P^/A~A]u;^^^ with a + (3 = 1 , (A3) 

which then leads to je Ao+'jo A^ = je+lo- This relation, together with the previous identity jgAo = 'YoAe, determines 
the values of the two amplitudes A^ and Ao- Demanding periodicity for the ansatz (A3) results in a second identity 
for the two integration constants a, /3, namely /3 = au^ . Comparing the result with our assertion, Eqn (A2), we 
identify ui = ui- and Z{uj) = 1 + ui^ . Since ui- <1, the N ^ oo limit is also obvious. 

To treat the general case, wo set up some notation. For any even m > 2, we define two auxiliary functions which 
depend on an ordered set of arguments 1 < fci < fc2 < ... < fcm < N, 

SiP{k,,k2,...K) ^ ^Ak, Afc, = Si^'\k,,k2)Si^\ks,k4)...Sif\km-l,km) 
Si^\ki,k2,...km) = VAk, Afe,...Afe^a;-('=^-'=^)-('=^-'=3)--('='"-'='"-) = Si''\kuk2)Si'\k3,ki)...Si'\km-l,km) 



Thus, we rewrite Eqn (A2) in the form 



Z{u;) 
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and recall that both S'^^\i,j) and S2'\i,j) satisfy the equations of motion for the infinite chain, i.e., (10), (11), and 
(15). 

We now proceed as follows: First, we argue that both Sm\ki,k2, ■■■km) and Sm{ki,k2, ■■■km) solve the homogeneous 

and inhomogeneous equations of motion (9) for m-point functions on the infinite chain. For Sm\ki,k2, ■■■km) = 

limTv^oo i'^kf-O'k^) , this is immediately obvious. For Sm {ki,k2, ■■■km), we simply need to retrace the inductive 

proof from Section 4: Since the basic relations (10), (11), and (15) all hold for \ we may replace each (afci...(7fe^) in 

Eqns (17-18) by the appropriate Sm without violating any of the equalities. Since some of the equations of motion (9) 

for the m-point correlations are inhomogeneous, only a convex combination aSm) {k\, ■■■km) ~t~ (^^rn {k\, ■■■km) 
with a + /? = 1 is a solution. 

Finally, we have to to satisfy the last remaining property of the full correlation function, namely, periodicity. Given 
an ordered set of arguments 1 < fci <k2 < ... < km < N,we assert that a general correlation function of m > 2 spins 
takes the form 



(cr/ji...cr/j^) = for odd m 



\ _ V^fci Afc2---^fcm f lk2-kl) + {k4-k3) + ... + {km-km-l) , , .AT- [(fes -fel ) + (fe4 -fe3) + . • . + (fem -fem- 1 )] \ 



'{OJ) 

^ {5({)(fci,fe,...fcn^) +t^'^^i'H^i,fc2,.../cm)} for cvcu m 



(A4) 



Z{uj) 



which is obviously periodic on the ring. Comparing Eqns (Al) and (A4), it is manifest how closely the correlations 
of the non-equilibrium model mirror those of the Ising chain. Again, the only significant difference is the appearance 
of the two amplitudes. A,, and Aq. 
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